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ABSTRACT: A rotor acoustic prediction computer program (RAPP) has been devel- 
oped to further the understanding of main rotor noise. The prediction program is based 
on a well known solution to the Ffowcs Willliams and Hawkings equation for acoustic 
pressure, wherein nonconr.pact monopole terms model rotor blade thickness and distributed 
dipoles model local blade surface pressure. RAPP has been validated by comparisons with 
the Boeing model 360 DNW/AATMR acoustic pressure base. This is the first work that uti- 
lizes both the acoustic pressure and blade surface pressures acquired during this test. Future 
applications of RAPP to the study of main rotor noise are discussed. 
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Introduction 


The rapid development of passive acoustic arrays and signal processing technology 
that detects, tracks, and identifies a rotorcraft poses a growing threat to the survivability of 
U. S. combat helicopters .✓Xn essential prerequisite to reducing the acoustic detectability 
of military rotorcraft, and hence to increase their survivabHitj^is a better understanding of 
main rotor noise,' which is the major contributor to the overall noise. A simple, yet accurate. 
Rotor Acoustic Prediction Program (RAPP) has been developed to advance the understand- 
ing of main rotor noise. This prediction program^tilize^the Ffowcs Williams ar.d Hawkings 
(FW-H) equation. The particular form of the FW-H equation* used in this analysi>is well 
suited for the coupling of the measured blade surface pressure to the prediction of acoustic 
pressure. 

The FW-H equation is an inhomogeneous wave equation that is valid in all space and 
governs acoustic pressure generated by thin moving bodies. The nonhomogeneous terms 
describe mass displacement due to surface motion and forces due to local surface stresses, 
such as viscous stress and pressure distribution on the surface. There is a rich history of the 
application of the FW-H equation as applied to the prediction of propeller and rotor noise. 
The events of the past, which have contributed to the wealth of knowledge in which this 
paper has its foundation, are discussed in depth in references 1 and 2. The history begins 
with Gutin’s 3 compact source assumption. He modeled steady aerodynamic forces on a 
propeller with harmonic acoustic dipoles in the frequency domain. Garrick and Watson 4 
extended this approach to uniform rectilinear motion. Denting, 5 Amoldi, 6 and Lyon 7 stud- 
ied the effects of blade thickness on acoustic noise respectively for static propellers, uniform 
rectilinear motion, and for helicopters in forward flight. Farassat, 1 Hawkings and Lowson, 8 
Isom, 9 and Schmitz and Yu 2 have used noncompact monopole terms to model rotor blade 
thickness and distributed dipoles to model local blade surface forces. r 

There aredour types of main rotor noise: BVI noise^ low-frequency noise* l)igh-speed 
impulsive noise, and broadband noise. 10 JBlade- vortex interaction noise occurs when a tip 
vortex, previously shed by a rotor blade, passes close enough to a rotor blade to cause large 
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variations in the Made surface pressures. This event is most disturbing when it happens 
on the advancing side of the rotor disk. ’^Low-frequency noise includes hover and low to 
moderate speed forward flight. For these flight conditions, the low frequency components 
of the acoustic signal dominate. RAPP was developed for BVI and low-frequency appli- 
cations* These are the flight regimes producing the propagating noise that leads to acoustic 
detection. ! ; : ■* j 

High-speed impulsive noise is a highly nonlinear, transonic phenomenon that is beyond 
the scope of this work. Recent efforts that do address high-speed impulsive noise couple 
computational fluid dynamics (CFD) with acoustics. The particular method of Isom and 
Purcell 12 couples nearfield pressure from a full potential finite-difference method with a 
Kirchhoff integral formulation to extend the CFD results to the far field. Broadband noise 
will not be addressed in this paper either. It is of comparatively small amplitude and is a 
result of unpredictable blade loading such as that caused by turbulence. 13 

RAPP implements the noncompact methods to efficiently and accurately predict BVI 
noise and low-frequency noise characteristics generated by the main rotor. This paper will 
explain the details of the blade thickness modeling and the blade surface force modeling 
used in RAPP. It will establish the validity of RAPP by comparisons with aeroacoustic 
experimental data for BVI and low-frequency noise. Future applications that are made pos- 
sible by prediction programs such as RAPP are discussed. The most important application 
is to use RAPP in conjunction with an aerodynamic blade surface pressure prediction pro- 
gram to leam to design quieter rotor blades. Hence, this would minimize detection ranges 
and decrease the threat to the survivability of combat helicopters. 

Development of the Rotor Acoustic Prediction Program 

As stated earlier, RAPP is based on the FW-H equation. This integral equation is 
expressed as follows: 
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Equation (1) is an integral equation for the acoustic pressure, p\ at a point x and a time 
t. The subscript ret indicates that the integrals are to be evaluated at the retarded time 
t - t — r/a o . The speed of propagation of sound in the fluid medium is designated by 
ao , and po is the density of the undisturbed fluid medium. M r is the Mach number in the 
radiation direction, M r = In this equation for the radiation Mach number, Vi is the 
local free stream velocity vector, and f, is the unit vector describing the radiation direction. 
For the above form of equation (1), the radiation Mach number should not equal 1. As a 
rule, the advancing tip Mach number should remain below the critical Mach number for the 
airfoil section. RAPP does not contain any nonlinear terms and should not be applied to 
nonlinear circumstances such as the transonic flow regimes. 



The first term on the light-htmd side of equation (1) is referred to as the thickness term. 
This term represents the disturbance of the fluid medium caused by the airfoil. For a fixed 
volume, a short fat airfoil will ca ase a greater displacement than a long slender airfoil. This 
is the influence of the normal ve’.c ; ty on the airfoil surface. Thin airfoil theory defines the 
normal velocity as 

t n = Vole* + tan" 1 (2) 

ax 


Here a is the airfoil angle of attack, Vo is the local freestream velocity, and g is the thick- 
ness distribution of the airfoil. Even though the angle of attack is not zero, the term 
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is negligible. It does not contribute to the thickness noise, even for a > 10°, and has been 
neglected by other researchers, namely Schmitz and Yu 2 . The term 
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can be approximated quite well by a chordwise distribution of sources and sinks at about 
five span wise locations. The sources arc assumed to be of strength 
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Equation (4) can be discretized as 
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Here dS = dxdz, and the subscript i; indicates quantities that arc functions of both the 
chordwise and span wise variables. The thickness of the airfoil at the ith spanwise location 
and the jth chordwise location is dyij. The thickness noise has been well modelled with 
10 chordwise sources at each of 5 spanwise radial stations. Adding more sources and more 
stations does not improve the accuracy of the thickness prediction. It just lakes more com- 
putational time. This was verified by adding more radial stations and more sources at these 
radial stations to RA PP in the dr elopmental stage. 

The second and third terms in equation (1) are the loading terms. Only the second 
term in equation (1) will be used to predict the loading noise for two reasons. The third 
term, which does not include a partial time derivative, is considered the near-field term. 
It is divided by the square of the distance between the source and the observer. For an 
observer in the far-field, this near-field term is negligible. Also, the quasi-steady near-field 



acoustic pressure was not measured by the microphones used in acquiring acoustic test data. 
These microphones were back vented free-held response microphones. Consequently the 
frequency of the quasi-steady near-field acoustic pressure was below the low frequency 
response of these microphones. 

RAPP uses the acoustic lifting line formulation to model the loading noise. If the 
loading information is supplied as blade surface pressure, the sectional lift and drag are 
determined by integrating the surface pressure over the chord. The loading term can be 
discretized as follows: 
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The sectional lift and drag are defined by l and d. The normal vectors h* and ry are defined 
parallel and perpendicular to the tip path plane o f "he rotor. The airfoil is exerting a force 
on the fluid, therefore, convention defines lift as positive downward and drag as positive 
toward the leading edge of the airfoil. Figure (1) displays the sign convention of these 
normal vectors in the blade fixed coordinate system. Two important assumptions that have 
been made in arriving at equation (7) are that the distance r and the radiation Mach number 
M r are held constant over the chord, c t . The loading terms are modelled well by equation (7) 
with five radial stations on the outer 50 percent of the rotor span. 



Figure 1 . Presentation of the sign convention for the surface normals in the blade Fixed coordinate system. 


The thickness and loading terms must be evaluated at the retarded time, r. This allows 
for the proper summation of the acoustic disturbances at the observer location and time. An 
iterative technique is needed to solve for the location of all the sources that contribute to the 
acoustic disturbance for a time t and observer x. Iteration is necessary to solve the re tardea 
time equation because the distance between the source and the observer is a function of 
the source time. The locations of the contributing sources comprise the acoustic planform. 
The acoustic planform is defined as the locus of emitting source locations whose signals 
arrive simultaneously at a time t to the observer. Please note that the surface over which 
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planform. 

Both the Newton iteration method and the Bisection method are used to solve for the 
source locations on the acoustic planform. Newton’s method converges quickly, or in just 
a few function equations, if the initial guess is “close enough” to the solution. Although 
the function on which the iterations are being performed, 


/(i) = i - (r- Ar+ t(t- Ar)/ao) (8) 

is a continuously increasing function, //( x ) can become very small. Since this term is the 
denominator in the iteration scheme, the method will not converge unless the initial guess is 
close enough to the solution. To guarantee convergence of Newton’s me..iod, a satisfactory 
initial guess is supplied by the Bisection method. This method is not affected by //( x ) , 
but it does require many function evaluations. This is why Newton’s method is used to 
conclude the iteration after the Bisection method finds a suitable first guess. 

A central difference method, which is second-order accurate, is used for the time 
derivative. The time derivatives of the thickness and loading terms are taken separately, 
so that the contribution to the acoustic disturbance from each phenomenon can be studied 
separately. 


Validation of the Prediction Code 


Discussion of Test Data 

RAPP was validated by comparisons with the Boeing 360 experimental acoustic test 
data acquired during the Aerodynamic and Acoustic Testing of Model Rotors (AATMR) 
program at the Duits-Nederlandse Windtunnel (DNW). The DNW Boeing 360 model rotor 
acoustic database has been reviewed and analyzed by Zinner et al., 10 and is considered 
of high quality to validate acoustic prediction codes. Zinner shows that the Boeing 360 
acoustic *"st data displays the expected trends for high-speed impulsive noise, blade-vortex 
interaction noise, low-frequency noise, and broadband noise. The conditions for maximum 
BVI noise were shown to range from low shaft-tilt angles at high advance ratios to high 
shaft-tilt angles at low advance ratios. BVI noise has maximum radiation in the direction 
toward the advancing side of the rotor and at least 25 degrees below the rotor plane. A 
severe test of RAPP’s predictive ability will be whether or not it can predict these trends for 
BVI. 

The accuracy of RAPP’s predictive ability depends on the accuracy of the input blade 
surface pressure. The DNW Boeing 360 experimental surface pressure data has been uti- 
lized by Caradonna et al* 4 in '.he prediction of airloads and is considered to be of the quality 
necessary to support test/theory. Thus, this DNW Boeing 360 blade surface pressure was 
thought to be a good source of input to validate RAPP. Both Caradonna and Zinner have 
discussed the DNW test facility, the data acquisition system, and the characteristics of the 
model rotor in great detail. The interested reader should study references 10 and 14 for 
more information on this experiment. 

The sectional lift and drag, supplied as input for RAPP, were obtained by integrating 
the Boeing 360 surface pressure over the chord. This technique did not work well for the 
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ficient number of pressure taps to get an accurate estimate of the sectional drag. Analytical 
estimates of the sectional drag are better suited for this application and are easy to deter- 
mine for hover where the loading is ideally constant about the azimuth. A good source of 
the sectional lift and drag is a potential flow lifting surface analysis program, such as that by 
Analytical Methods Incorporated (AMI). 15 The sectional lift and drag from AMI’s hover 
code was used in RAPP to predict the helicopter noise in hover. For the forward and de- 
scending flight regimes, the DNW Boeing experimental surface pressure data was used as 
input. The technique of integrating the surface pressure chordwise to get the sectional lift 
and drag worked well for these flight regimes because the lift dominates the loading noise 
contribution. 

Figure 2 shows the microphone areay that was used to acquire test data. These mi- 
crophone locations define the observer positions. This array of microphones is located a 
distance of three rotor radii from the rotor hub. 



Figure 2. Microphone array used 10 acquire acoustic pressure for the Boeing model 360 DNW/AATMR test 



RAPP was first applied to predict acoustic pressure in the hover flight regime. Since 
the loading is ideally constant around the azimuth for the hcvering rotor, it was a natural first 
step in the developmental process. Five hover test points wtre chosen from the Poeing 360 
data base for the hover comparisons. All five test points had the same solidity weighted 
thrust coefficient Ct/o = 07 , with hover tip Mach number ranging from Mh = .5 to 
Mh = 636 . The loading, or the sectional lift and drag, from the AMI Hover code supplied 
input to RAPP to match these five test conditions. In Figure 3a, the plots of predicted 
peak negative pressure versus Mach number depict the trend for hover of decreasing peak 
negative pressure below the rotor lip path plane. The peak negative pressure is defined as 
the most negative peak amplitude in the acoustic pressure time historic' At this loading 
condition, the thickness noise is dominant and its amplitude decreases below the tip path 
plane. The loading noise, which becomes predominant below the plane, does not increase 
more than the thickness decreases. Predicted peak negative pressure versus Mach number 
plots are compared to the measured plots of peak negative pressure versus Mach number 
in figures 3b, 3c, and 3d. Figure 3b shows comparisons for the in-plane microphone, d> = 
0°. The comparisons for the in-plane microphone are quite acceptable. Figure 3c shows 
comparisons for the microphone at $ = 15° below the tip path plane, and figure 3d shows 
comparisons for the microphone at 0 = 25° below the tip path plane. The comparisons in 
figures 3c and 3d are not acceptable. The experimental plots do not display the expected 
trend of decreasing below the tip path plane. This is due to recirculation of the air in the 
DNW test chamber. The DNW is not a hover chamber and does not allow for clean hover 
test conditions. A mean flow and recirculation occur in the test chamber and cause unsteady 
effects in the loading and in the acoustics. These unsteady affects are most noticeable in 
the out-of-plane observer locations and are quite noticeable in the plots shown in figures 3c 
and 3d. However, the in-plane comparisons, as shown in figure 3a, are good enough to 
conclude that RAPP can accurately predict acoustic pressure for the hover flight regime. 

The next step in the development of RAPP was to predict acoustic pressure in forward 
flight and in descent where BVI occurs. In these flight regimes, RAPP will have to contend 
with blade loading that varies about the azimuth. To get an accurate prediction, the blade 
loading must be known at many locations around the azimuth. The measured experimental 
pressure is known at 1024 locations around the azimuth. This resolution of almost every 
one-third of a degree enables the acoustic prediction program to compute accurate pressure 
time histories. Computer programs that calculate blade surface pressure for forward flight 
and BVI flight regimes are beginning to have the ability to predict the blade loading at very 
fine increments in azimuthal angle. The application of this capability to the Boeing 360 
rotor will be discussed in the application section of this paper. 

To validate RAi . or the forward and descending flight regimes, the measured surface 
pressure data from the Boeing 360 mode! rotor test was used as input. The surface pressure 
was integrated over the chord to get the sectional lift and drag at each radial station and az- 
imuth angle. The acoustic pressure time histories, predicted using the experimental surface 
pressures as input to RAPP, compare well with the measured acoustic pressure time histo- 
ries. This is presented in figures 4 and 5 which show comparisons of predicted to measured 
acoustic pressure time histories. Figure 4a shows the in-plane comparisons for a moderately 
high speed forward flight test point. Figure 4b shows comparisons for observers 25° below 
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the tip path plane for this same test point Figure 5a shows the in-plane comparisons for a 
descent or BVI test point. Figure 5b shows comparisons for observers 25° below the rotor 
plane. This data for BVI portrays the expected trends where the highest peak-to-peak levels 
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Figure 5b. Plots of acousuc pressure umc histones for microphones 25° below the tip 
path plane. The flight conditions are |i = 0.225, CT/ct - 0.07, and shaft tilt = 4.38°. 



are near 150° azimuth and 25° below the plane. It can be concluded from the comparison 
plots in figures 4 and 5 that RAPP is capable of predicting accurate acoustic pressure time 
histories for moderate-speed forward flight and BVI. 

Applications 

A prediction program such as RAPP, can be applied to study many aspects of main 
rotor noise. The contribution of the three noise generating mechanisms — thickness, lift, 
and drag — to the total acoustic pressure disturbance can be plotted separately. This allows 
for a more in-depth study of each noise generating mechanism. An example of such a plot is 
shown in figure 6 for an in-plane microphone from the same forward flight test point shown 
in figure 4a. 

The prediction of acoustic pressure is dependent upon the knowledge of the blade 
loading. Reliable sources of predicted blade loading are becoming available as the research 
in blade surface pressure prediction progresses. The most promising technique developing 
is the coupling of a full-potential rotor flow solver with a comprehensive lifting surface 
analysis. If is becoming possible for the techniques to produce the surface pressure at the 
number of azimuthal angles necessary for accurate acoustic prediction for forward flight 
and BVI. A very fine resolution in azimuth angle is necessary for accurate BVI predictions 
because the surface pressure changes so abruptly in a short source time increment Possible 
future applications for RAPP are to use the blade surface prediction capabilities to study 
the effects of varying airfoil sections, twist distributions and planform shape on BVI. This 
application could lead to helicopter rotor systems that are designed to be quieter in the BVI 
flight regime. 

Another application of RAPP, when used in conjunction with blade surface pressure 
prediction codes, is to predict experimental wind tunnel test data before the test is run. This 
would allow the experimentalist to decide the best microphone specifications and locations 
and at which test conditions to run. These predictions of the test data could help the exper- 
imentalist decide on the integrity of the data being acquired. The overall effect would be 
better planned wind tunnel tests that would provide more useful data. 

Conclusions 

The most important conclusion that can be drawn from the information presented in 
this paper is that RAPP can predict accurate rotor acoustic pressure time histories. The no- 
table technique and assumption used in RAPP is the acoustic lifting line formun tion. This 
formulation works well and considerably simplifies the application of the FW-ri equation. 
RAPP is a simple enough program to run on a personal computer and executes in RAP- 
Pidly. Therefore, the above mentioned applications could take place in the test field or in 
an industrial setting where supercomputers are not the norm. 

Some interesting byproducts of RAPP’s developmental process concern the details 
of the computer program. Details such as 10 sources along the chord at 5 radial stations 
being necessary and sufficient to accurately predict thickness noise. It is also interesting 
that even BVI noise can be predicted accurately by modeling the loading terms at only 
5 radial stations. This indicates that the azimuthal resolution is more important than the 








spanwise resolution. The fact that the acoustic predictions for hover rely so heavily on 
accurate section drag input will prove useful in future hover acoustic predictions. 
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